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I. INTRODUCTION 



The gauge invariant field strength vacuum correlators 

{Ff jLlUl (xi)U(x 1 ,X2)F^ 2U2 (x 2 )U(x2,x z ) ■ • - F llnVn (x n )U(x n ,xi)), (1.1) 

rx 

where U(x,y) = Pexp ig / dz >1 A fM (z) is the Schwinger color string, play a relevant role 

Jy 

in gluodynamics with and without quark sources. We know that in the infrared region 
these correlators are dominated by their non-perturbative behaviour. In particular the non- 
perturbative "gluon condensate" 

cy 

(^F« u (0)Fr(0)) noRpcit = F 2 , (1.2) 

7T 

plays a crucial role in the QCD sum rule method [p. 

The non-perturbative part of the gauge invariant two-point field strength correlator 
(F flu (x)U(x,0)F p(T (0)U(0,x)) has been calculated on the lattice, with the cooling method 
in and in presence of sources in ||. We define (in Euclidean space-time, as in the rest of 
this work) the gauge invariant correlator H 



(g 2 Fp U (x)U(x, 0)F Ap (0M0, x)) = (<V<W - 6 fip 5 I/X )g 2 D(x 2 ) + l - 



{x\8 up - x p 8 v \) 



dxp 



+ 7^- (x p 5px - x x Spp) 



g 2 D 1 (x 2 ). (1.3) 



A parameterization of the form 



D{x 2 ) = Ae~^ + - d e~^ D 1 (x 2 ) = Be~^ + - d e~^ T " (1.4) 



x 4 v ' x 4 



A ~ 128 GeV 4 , B ~ 27 GeV 4 , a ~ 0.69, b ~ 0.46, 
T g ~ 0.22 fm, T p ~ 0.42 fm, 

yields a very good fit to the (cooled) lattice data @ in the range 0.1 fm < x < lfm. At 
short distances the 1/x 4 term, which is of perturbative origin, dominates, while at distances 
x > 0.4 fm the non-perturbative term, proportional to e~ x l Tg , becomes more important. In 
an Abelian theory without monopoles the Bianchi identities yields D — 0. 



In the Stochastic Vacuum Model (SVM) j4|-[F] it is assumed that for processes which 
can be reduced to the calculation of Wilson loops with quasi-static sources (such as heavy- 
quark potentials and soft high energy scattering amplitudes in the eikonal approximation) 
the infrared behaviour of QCD can be approximated by a Gaussian stochastic process in the 
field strength and is thus determined approximately by the correlator (|1.3|) . In particular 



also the Wilson loop average is given only in terms of (1.3). 

We know from the strong coupling expansion and lattice simulations that the Wilson loop 
is an order parameter for confinement. The confining area law behaviour of the Wilson loop 
is reproduced by the stochastic vacuum model provided that the form factor D is different 
from zero and is dominated in the infrared region by a decreasing behaviour with the fall off 
controlled by a finite correlation length T g . These features of D are compatible with lattice 
data (see ( |1.4j) ). Furthermore this model gives a good description of certain features of high 
energy scattering (e. g. ||). We will come back to this point in Sec. 2. 

It is the goal of this paper to relate the gluon correlator to the Mandelstam-'t Hooft 
dual superconductor mechanism of confinement |J. In this picture the physical essence of 
the confinement is the formation of color-electric flux tubes between quarks due to a dual 
Meissner effect. The monopoles condense and lead to a dual superconductor which forces the 
color-electric field lines in flux tubes which are the dual analogue to the Abrikosov-Olesen 
strings. The formation of an electric flux tube is also the consequence of the Stochastic 



Vacuum Model 1C 



Furthermore, in an Abelian projection of QCD, monopoles are the degrees of freedom 
responsible for confinement. Monopoles condensation has been observed on the lattice (for 
a review see |ll|) and when confinement can be derived analytically (compact electrody- 
namics, Georgi-Glashow model, and some supersymmetric Yang-Mills theories), it is due 
to the condensation of monopoles. The monopole potential can be measured in the Abelian 



projection and it turns out that in the confining phase it has the Higgs form ||12|| . Lat- 
tice measurements of the distribution of monopole currents indicate that at large distances 
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gluodynamics is equivalent to a dual Abelian Higgs model, the Higgs particles are Abelian 
monopoles and these are condensed in the confining phase. In the maximal Abelian gauge 
the structure of the interquark flux tube was intensively studied and high precision measure- 
ments of the colour fields and the monopole currents recently allowed for a detailed check of 
the dual superconductor scenario with respect to the Ginzburg-Landau equations [13]. 



Analytic models of the infrared dynamics of dual QCD with monopoles were constructed 
[ |l4]Jl5| and their phenomenological consequences intensively investigated. In the effective 
dual model of Baker, Ball and Zachariasen (dual QCD, DQCD), the complete semirelativistic 
quark- ant iquark potential, the flux tube distribution and the energy density were obtained 
from the numerical solution of the coupled non-linear equations of motion and compared 
very favourably with recent lattice data |16|-p^]. Although the Lagrangian of this effective 
dual theory for long distance QCD is based on a non-Abelian gauge group, the results for 
the qq potentials aside from an overall color factor can to a very good approximation be 
described by a (dual) Abelian Higgs model. Therefore, the results are in this case largely 
insensitive to the details of the dual gauge group and the quarks select out only Abelian 
configurations of the dual potential RI5f| . 

Since an effective Abelian description of the infrared confining dynamics of QCD (at 
least for heavy quarks) emerges either from QCD (via Gaussian approximation and bilocal 
strength tensor correlators^) or via an effective Abelian Higgs models it becomes extremely 
interesting to exploit in which sense the two Abelian descriptions are equivalent and, once we 
assume an equivalence, what kind of constraints this imposes on the form of the QCD field 
strength correlators. In the present work we will obtain from the dual Abelian Higgs model 
information on the form of the gauge invariant two-point field strength correlator ( |1.3|) and 
in addition we will obtain an analytic approximation for the static heavy quark potential 



1 In the treatment of two Wilson loops, however, the non-Abelian characteristics of QCD becomes 
very important (||,|l0| . 
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given by the dual theory. 

There are two more arguments which motivate such an investigation. First, as we will 
discuss briefly in Sec. 2, the existence of a non- vanishing form factor D in the two-point field 
strength correlator of QCD seems to suggest quite naturally the existence of an effective free 
dual Abelian theory "behind" the long-range dynamics of QCD. Second, a recent comparison 
between the complete semirelativistic potentials obtained in DQCD and in the Gaussian 
stochastic approximation of QCD in the limit of large interquark distances showed up quite 
striking and surprising similarities (see @). The following analysis wants to shed some light 
on that. 

The plan of the paper is the following. In Sec. 2 we recollect some essential features of 
the gauge invariant two-point correlator in QCD and we establish a relation with the Wilson 
loop. In Sec. 3 we investigate the analogous quantity in the dual Abelian Higgs model 
without sources. For a constant Higgs field we reproduce a two-point correlator having the 
same behaviour as obtained by other authors by studying the London limit of a dual Abelian 
Higgs model. In Sec. 4 we introduce sources and obtain an analytic expression for the static 
potential. This suggests a connection between the parameters of the two-point field strength 
correlator in QCD and those of the dual Abelian Higgs model. Finally, Sec. 5 contains some 
conclusions. 

II. GAUGE-INVARIANT TWO-POINT GLUON CORRELATOR AND WILSON 

LOOP IN QCD 

We consider the correlator of two gluon field strengths in QCD at different space-time 
points, connected by a Schwinger string. This string can either consist of two strings in the 
fundamental representation or one string in the adjoint one. For definiteness in notation we 
choose the first possibility and consider the quantity 

(g 2 F^(x)U(x,0)F Xp (0)U(0,x)). 
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The Lorentz decomposition of this correlator is given by Eq. ( |1.3| ) and the results of the 
lattice measurements are collected in Eq. ( |1.4| ). The leading (tree level) perturbative contri- 
bution is contained in the form factor D\. In an Abelian gauge theory without monopoles 
the Bianchi identity implies that the form factor D vanishes identically M. In a non- Abelian 
theory or in an Abelian theory with monopoles D can be different from zero. Let us briefly 
review how a non- vanishing D leads to confinement 

In the presence of heavy quark sources the relevant object in QCD is the Wilson loop 
average H^(T), where T is a closed curve built up by the trajectories of external sources 
and some Schwinger strings connecting the end-points. By means of the non-Abelian Stokes 
theorem ]2TJ] one can express the Wilson loop average H^(r) in terms of an integral over a 
surface S(T) enclosed by the contour T. A way to evaluate analytically this quantity consists 
in expanding this expression via a cluster expansion and keeping only the bilocal cluster (i. 
e. in assuming that the vacuum fluctuations are of a Gaussian type)[|: 

W(T) = exp (ig jf dz,A^ } ^ ^P exp (ig jf ^ dS^F^u, s )) \ (2.1) 

Q ^exp(-i/ dS pu (u) I dS Xp {v)(g 2 F flu (u,x )F Xp (v,x ))) , (2.2) 
svm y 2 Js(r) Js(V) J 

where P F^u, xq) = P U(u, xo)F fJ/V (u)U(xo, u). Assumption (|2.2|) corresponds to the so- 
called Stochastic Vacuum Model (SVM) [Q. The point x Q is an arbitrary reference point on 
the surface S(T) needed for surface ordering. Of course the final result in the full theory does 
not depend on the reference point xq. The results obtained in the Gaussian approximation, 
however, will generally depend on it. This dependence is minimized by choosing S(T) to 
be the minimal area surface with contour T ||22|| . Under this condition one may neglect the 
Xq dependence on (g 2 F pjV (u,x )Fx p (v,x )) and recover in this way translational invariance. 

2 For an extensive discussion on the validity of this assumption see j/J. Moreover, recent lattice 
calculations seem to confirm that heavy quark potentials are really dominated by the two-point 
gluon field strength correlator Q. 
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Then, the decomposition of Eq. ( |1.3| ) can be used (by replacing x 2 with {u — v) 2 ). 

All the spin and velocity dependent potentials up to order 1/m 2 in the quark mass can 
be expressed in terms of the functions D and Di P, |2"0|J2"B[ • In particular the static potential 
is given by 

V (R) = £ / d 2 x (R - \ Xl \) D{x 2 ) + ^D^x 2 ), (2.3) 

Z J\ Xl \<R Z 

with d 2 x = dxidx^, x 2 = x\ + x\. The string tension emerges for large qq distances R as 

2 

a = !L J d 2 xD(x 2 ). (2.4) 



Therefore a no n- vanishing D function leads to confinement .0 

While lattice data confirm the existence of a non- vanishing form factor D with exponential 
fall off, up to now there is no analytic tool which allows to calculate and to interpret the 
non-perturbative contributions to D in the long-range regime. 

We observe, however, that a non-vanishing function D emerges naturally if we assume 
that there exists an effective "dual" Lagrangian describing an Abelian gauge theory for which 
the dual two-point field strength correlator coincides in the long-range limit with the bilocal 
cumulant given by Eq. (|1.3|) . Let us call Gp V the (Abelian) field strength of the dual theory. 
Since we assume this theory to observe the Bianchi identities we have in general 



(g 2 G, u (x)G Xp (0)) = \ 



d d 

(x\5 vp - x p 5 u x) + - — (x p 5px - xxSpp) 



OX a P ' OX, 



9 2 A(x 2 



The expectation value of the dual of the dual fields, Gp U = -tp^pGap, is 

2 



3 The 1/x 4 term in D in (1.4) is a one-loop perturbative contribution [24] and has not to be consid- 
ered in the calculation of the string tension a. Preliminary results indicate that these perturbative 
contributions to D appearing at one loop and higher orders are cancelled by higher order correlator 



contributions 2E]. This is not surprising since in a non- Abelian theory perturbative contributions 
beyond the tree level are surely not of a Gaussian type. 
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(g 2 G pi/ (x)G Xp (0)} = (d&Sup - 5, p 5 uX )g 2 (a(x 2 ) - x 2 -^A(x 



+ 2 



dx 

d d 
— (x\6 vp - x p 8 u x) + - — (x p 5 mA - xx&w) 



g 2 d d 2 
dx 2 dx 2 



It shows a tensor structure like the one multiplying D in equation ( |1.3| ). The existence of 
such a correlator therefore seems to suggest the existence of a dual Abelian gauge theory for 
which at big distances the field strength correlator behaves as the corresponding correlator 
of the dual theory: 

(g 2 F pu (x,x )F Xp (0,x )) ~ (g 2 G^(x)G Xp (0)). (2.5) 

In the next section we want to explore some consequences of Eq. fl2.5| ). In Sec. 4 Eq. 
( ^.5| ) will be replaced by a better founded assumption on the Wilson loop. Nevertheless the 
basic idea behind (|2.5| ) will remain. 



III. DUAL ABELIAN HIGGS MODEL WITHOUT SOURCES OR VORTICES 

The aim of this section is essentially pedagogical. We will reproduce in a clear and 
economical way some of the results which can be found in the existing literature on the 
London limit of a dual Abelian Higgs model. We will prove in this way that assumption 
( P-5|) is reasonable, i. e. compatible with ( |1.4|) . We will also show the drawbacks of this 
approach and try to justify why we need to take into account external charge sources. This 
will lead to the results of Sec. 4. 

Let us consider a very naive context, i. e. a "dual" vector gauge field C p minimally 
coupled with some external scalar field which we could call a Higgs field. The action is 
given by 



S{C P ,4>)= d 4 x 



\G pu (x)G pv (x) + i( J D M 0)*(x)( J D M 0)(x) + V{<f>*(x)<f>(x)) 



(3.1) 



where G pv (x) = d p C v (x) — d u C p (x) and V(</>*0) = -^(4>*4> — 4>l) 2 (with 4>o different from 
zero). The Higgs field is coupled to the gauge field C p via the covariant derivative D^cj) = 
(d p + ieCft)<p. 



We choose a gauge in which the regular part of the phase of vanishes, the so-called 
unitary gauge. The propagator, JCp U = (CpC u ), of the field C p satisfies the equation: 

(d 2 5 u , - d u d, - e 2 <f ) 2 (x)5 u ,))C ua (x, y) = -5, a 5\x - y). (3.2) 

The quantity in which we are interested is what we could call the "dual" of the field strength 
two-point correlator in the theory described by the action (|3.1| ): 



Gay\ P (x, y) = (5\ a 6fry - 5 x -yS pa )5 4 (x - y) - e^Xpepa^d^JCu^x, y). (3.3) 

For a matter of convenience we prefer to define Q a ^\ P with the delta contribution subtracted 
out explicitly. In this model Q al \ p is the equivalent of the quantity (g 2 G flu (x)G\ p (y)) intro- 
duced at the end of the last Section. Eqs. ( |3.3| ) and (|2.5|) then give the correlator ( |1.3j ) in 
terms of the propagator of the dual theory. 

Let us study, now, the case where the Higgs field has the constant value 0o- Then, Eq. 
( fj.2| ) can be written as: 

(d 2 ~ e 2 4)}C™(x, y) = - (s^ - ^ S\x - y). (3.4) 

This is simply the equation defining the free propagator of a massive vector boson with mass 
M = e0 o : 



with/C°°(x-y) = j 



d 4 p „_ ip(x _ y) 1 M K x {Mx) 



-e 



(2tt) 4 p 2 + M 2 (2tt) 2 x ' 

where K n (n = 0, 1, 2, . . .) are Bessel functions. As a consequence we can write 



Ga~y\ P (x - y) = (<5xJ py - Sx 1 S pa )D°°((x - y) 2 ) + - 



d 

((x - y) a 5 pl - (x - y)~ l 8 P o) 



d 

+ -q^- \\x - yjjSxa - (x - y)Jx 1 ) 



dx x 

D?((x-y) 2 ), (3.5) 



with 
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D°°(x 2 ) = 5\x) + d 2 lC^(x 2 ) = M 2 JC°°(x 2 ) = Ml K ^ Mx \ (3. 6 ) 

An 2 x 



d ^^) = -aA-ic oo (x 2 ) }l 



dx 2 x 2n 2 x 2 



K ^ Mx K^MMx) + K 2 (Mx)) 



(3.7) 



x 2 

Therefore, the assumption that Q al \ p has the same long-range behaviour of the gauge in- 
variant two-point field strength correlator in QCD (see Eq. ( |2.5|) ) is compatible with the 
parameterization and leads to a correlation length T g equal to the inverse of the dual 
gluon mass M. In particular, the asymptotic behaviours of /C°° are: 

1 1 

1*1—0 (2ir) 2 x 2 



^(x 2 ),-^ -—;- + ■■■, (3.1 



K,°°(x 2 )^- f * -i , e- Mx + ---. (3.9) 
V (2tt) 3 / 2 ^/Fx 3 / 2 V ; 

The results shown here coincide with those obtained from the London limit of a dual 



Abelian Higgs model in p6[ . The seeming difference for what concerns the function D is 
due to the fact that we have subtracted out explicitly in our definition of Q al \ p the delta 
singularity which in the referred work is taken into account in a regularized form. One may 
wonder how the result of a topologically trivial model (no singular Higgs phase) agrees with 
results which take into account properly the internal Abrikosov-Nielsen-Olesen strings. This 
is due to the fact that G a -y\ P is sensitive to the string only via the strength of the Higgs field 
and this is fixed to a constant here as well as in the London limit of a dual Abelian Higgs 
model. 

The agreement between both approaches reveals a common weakness: the missing treat- 
ment of the interaction between the internal strings present in the dual Abelian Higgs model 
and the string between external quark sources. In Ref. f26[ no external sources were intro- 



duced and the result (|3.6| )- (|3.71 ) for the correlator was obtained in the following way. The 



functional integral for the Abelian Higgs model was rewritten in such a form as to exhibit 
integration over the closed surfaces of the (internal) strings. From the form of the contribu- 
tion of a single closed surface in the London limit it was deduced that it could be obtained 
by a correlator like (|3.6| )-( [OD in the Gaussian approximation. Implicitly this form was as- 
sumed to be valid also for external sources. When external quark sources are introduced 
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however, the strings of those sources will interfere with the internal strings. Some aspects 
of that phenomenon (the linking number) have been treated in P7| , but to our knowledge 



there exists no analytic attempt to evaluate the influence on the phenomenologically relevant 
parameters. 

We notice that, due to the short-range behaviour ( |3.8|) , Eq. Q3.7D reproduces the expected 
short-range behaviour of the function D\ (~ 1/x 4 , see ( p..4|) ). Due to the short distance 



behaviour of the function D, like 1/x 2 , the string tension we obtained using (2.5), Eq. (|3] 



and Eq. (^4), is logarithmically divergent: 

a oo _ 9\ f £x poo^ = w(j) 2 l imKo ( e ) „ ^2 (lQg2 _ log£ _ 1Q) 

where we have used the Dirac quantization condition e = 27r/g, relating g to the coupling 
constant e of the dual theory. The divergence is a short- distance effect and appears to be 
a result of the freezing of the Higgs field to the vacuum value <po, i. e., in terms of the 
dual Abelian Higgs model, of the London limit. Assuming a coordinate dependent Higgs 
mass going to zero like \x\ near the origin, would yield a finite short-range behaviour of the 
function D°° while preserving the perturbative short-range behaviour of the function Df. 
There is, however, no motivation for such an anisotropic behaviour of the Higgs field unless 
we introduce some charges into the vacuum. Only in such a context we can expect that near 
the sources and on the connecting flux tube string the Higgs field vanishes while far away it 
assumes the vacuum value 4>o- This will be precisely the subject of the next section, where 
we will consider a dual Abelian Higgs model with external charges and where we will also 



change our intuitive duality assumption (|2.5| ) to a more physically justified one. Moreover, 



we recall here that recent lattice data |T3| confirm that in the presence of external quark 



sources the distribution of electric fields and monopoles currents does not fulfill the London 
limit. 

To conclude this section we comment briefly on the translational invariance of the con- 
sidered correlators. As long as (ft is considered as an external field in Eq. ( |3.2| ), Q al \ p is not 
translational invariant and therefore in order to take advantage of the decomposition (|1.3|) we 
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have to fix our reference frame in such a way that the point y coincides with the origin. This 
fact is by itself not in contradiction with the duality assumption ( |2.5| ) since also the correla- 
tor in the direct theory, (F pu (x,xo)F\ p (y,x )), is in general not translational invariant, and 
only by choosing the reference point Xq on the straight line connecting x with y is invariance 
recovered. Finally, we notice that Gcry\p is translational invariant in some particular cases: 
if we assume <fi constant, as we have done in this section, or partially (in the longitudinal 
coordinates) if we assume that <fi depends only on some (transverse) coordinates. This last 
situation will be exploited in the next section. 



IV. DUAL ABELIAN HIGGS MODEL WITH EXTERNAL QUARK SOURCES 

In this section, for the reasons stated above, we want to consider a dual Abelian Higgs 
model with external quark sources. In particular we want to make a duality assumption on 
the long-range behaviour of the Wilson loop associated with the dynamics of a two heavy 
quark bound state. This assumption will take the place of our previous statement ( |2.5| ). We 
will see that some general features will, nevertheless, be preserved. 

Following we assume that the long-range behaviour of the Wilson loop average H^T) 
associated with a two heavy quark bound state is described by the functional generator of a 
dual Abelian Higgs model with external quark sources : 

W(T) ~ (e- 5(c -^), (4.1) 

where the bracket ( ) means the average over the gauge fields C p and the Higgs field <fi. The 
Abelian Higgs model is dual in the sense that it is weakly coupled. Therefore the right-hand 



side of Eq. (3~T) can be evaluated via a classical expansion. 

The action 5" is given by equation ( |3.1| ), but since we want that in the A = limit S 
describes the dual of a U(l) Yang-Mills theory with two external point-like charge sources 
— g (particle) and g (antiparticle), we define the field strength tensor G^, now, in such a way 
that it contains not only the dual gauge fields C M but also the field of the external sources 
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2§: 



G^{x) = d^C u {x) - d u C^ (x) + G s (x) 



(4.2) 



where 



rl rl d]Jc£ dy 'ft 

G lu( x ) = 9^^ al 3 J o dr da £ 5 (x-y(r,a)), 



(4.3) 



and ?/ m (t, a) is a parameterization of a surface S(T) swept by the Dirac string connecting the 
charges — g and g. Therefore S(T) is a surface with a fixed contour given by T (y^r, 1) = zi^ 
and y M (r, 0) = z 2 ^, where Zi^ and z 2ai are the charge source trajectories). Notice that the 
divergence of the dual of G is just the current carried by a charge g moving along the 
path V: dpG^p(x) = dz a 5 4 (x — z). The charge g is related to e by the usual Dirac 

quantization condition e = 2it/g. 

The leading long distance approximation to the dual theory is the classical approximation 



-s(c,4>) 



) ~ e 



-S{C$,4> oX ) 



(4.4) 



where C!f and <fi cl are solutions of the equations of motion: 



(d 2 5 Ufl - d u d^ - e 2 4> 2 (x)S ulx )C u (x) = -d v G%{x), 

(<9 M + ieC M (x))(<9 M + ieC^x))^) = A(0 2 (x) - <fo)<t>{x) 



(4.5) 
(4.6) 



Using these equations it is possible to write S{G 



as: 



+ J d 4 x 



-(d0(x)) 2 + \/(0 2 (x)) 



(4.7) 



where the propagator K va was defined by Eq. (|3.2|). Finally, integrating by parts, we obtain 

,2 , , rl 



-Xd(t){x)Y + V{cf)\x)) 
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where the tensor Q al \ p is the same as given by Eq. (|3~3|). Comparing with Eq. ( |2.2| ) we 
conclude that Q al \ p plays the same role in the dual theory as the two point correlator in 
the Stochastic Vacuum Model if we neglect the contribution of the Higgs field to the action 
in Q4.8| ) . In the London limit the contribution of the Higgs field to Q4.8] ) vanishes and the 
identification is exact. 

In the general case we are considering here also the Higgs part gives a contribution to 
the non-perturbative dynamics. But let us neglect the dependence of the Higgs field, via the 
equations of motion, on the strings and take into account the contribution coming from the 
Higgs part as a finite contribution to the string tension. Then, also in the general case, Q al \ p 
can be considered equivalent to the QCD two-point non-local condensate and in principle 
gives information on the validity of the decomposition ( |1.3|) and on the existence and the 
behaviour of the D and D\ functions. 

Notice that in the derivation of Eq. ( ^4.8| ) we have not considered surface-like contributions 
which would arise from the functional integral on the right-hand side of Eq. ( |4.1| ) once 
singular Higgs phase contributions are taken into account (these are also called Abrikosov- 
Nielsen-Olesen strings). These surface terms would interfere with the surface terms coming 
from the external quarks loop. We make the assumption that these interference terms are 
unimportant in order to evaluate the long-range behaviour of the (heavy quark) Wilson loop 
average after the duality assumption ( |4.1|) . In this way all the contributions coming from 
the singular Higgs phase factorize in the functional integral to a constant and play no role in 
the dynamics (see also the discussion on this assumption made in the context of the London 
limit in Sec. 3). 

We now evaluate ( |4.8|) beyond the London limit. 

Let us write a point x in the four dimensional Minkowski space as x = (x\\,x±), where 
) and x± = (x2, £3) are now two-dimensional vectors. Let us indicate with small 
letters the components of x belonging to x» (e. g. x a ,Xb, ■■■) and with capital letters the 
components of x belonging to Xj_ (e. g. xa,%b, •••)• 111 order to simplify the problem and to 
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allow us to give an analytic evaluation of (|4.8|) we choose the surface S(T) (see Eq. (fOD) to 
belong to the plane x± = 0. It is reasonable in this case to assume, at least far away from 
the charge sources (i. e. in the middle of the flux tube), that the Higgs field depends only 
on the transverse coordinate x±: 

= 0(x ± ). (4.9) 

We will make this crucial assumption for the rest of this section. 

From Eqs. ([OP and flOP we have K m (x,y) = JC^ a (x\\ - y\\,x±,y±). 
In this situation we have that Eq. ( |4.8|) can be written as: 

2 

S(C$, cl ) = S - I dS lA {x\\) [ dSu{y\\)Guu{x\\ - y\\) + Higgs sector, (4.10) 
Guu(x\\ ~ y\\) = S 4 {x\\ - y\\) - euAB€ucDd y c &X K, B d{x\\ - V\\,x±,y±) _ . (4.11) 

x_\_ — y_\_ — u 

After some simple manipulations it is possible to obtain from Eq. (|3.2|) an equation only for 
the transverse components of the gauge field propagator: 

d 2 ± 5 C B - d B d c - e 2 (j) 2 (x ± )5 C B\ K>ca{x\\ ~ V\\,x±,y±) 
+ d\ [5 CB - d B (d 2 x - e^^x^dc] JC CA (x\\ - y h x ± ,y ± ) = -S BA S\x - y), (4.12) 

where <9jj = d a d a and d\ = 8a9a- 

We look for a solution of Eq. (|4.12| ) of the type: 



eiACDdc JC B d{x\\ - y\\,x±,y±) n = -eucBXcJC(x\\-yn,x±). (4.13) 

This is reasonable since in the transverse plane we have rotational invariance. The function 
/C is unknown, but from Eq. ( [4.12j ) we have that it satisfies the equation: 



(<9 2 - e 2 2 (x ± ))x A /C(x||,x ± ) = -5 2 (x {l )d A 5 2 (x ± ). (4.14) 

In the limit for x± — > 0, we look for a solution xa^C of the type: 

x A JC(x\\,x±) = d A K?(x) + x A f(x\\)g(x±), (4.15) 
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where /C p is defined by 



therefore KJ 



d 2 K? = -5\x) 
1 1 



(2tt) 2 x 2 

and we — e , by mp0smg - , T h e unknOWn fumfas , md , satlsfy 

the equation: 

d\\fi x \\) (x A g(x±)) + f(x\\)dl (x A g(x±)) - e 2 2 (x_L)/(x||) (x A g(x ± )) = 

e 2 <f) 2 {x ± )d A IC p {x). (4.16) 

Integrating over the longitudinal coordinates both side of the equation, we get 

(d]_ - e 2 4> 2 (x i _)) (x A g(x±)) = -^-e 2 4> 2 (x ± )-^-, 
v 7 2n x L 

d 2 x\\d A K?(x) = 5-. This is exactly Eq. (JA5I) of the Appendix. 

2ir x ± — 

Moreover, also the boundary conditions are the same, since 

c n( x ) = J d 4 yK m (x,y)d v Gl a (y). 
Therefore, a solution exists (for small x±) and is given by 

g(x±) = - (4.17) 

Z7T X± 

For the definition of C np see the Appendix. Using the expansion ( JA7|) , for small x± we have: 

= — - — x ± + ■ • • , 

Z7T 1D7T 

= 7T~X A H , 

C2 

^XAg(xi) = --^ x A H , 

where and S 1 ^ are some constants defined as: S c = ]im^eC Dp (x ±) / x ± and = 



lim e<fi(x±)/x±. By solving numerically the static equations of motions (|4.5| ) and ( |4.6| ) 

xj_ — >0 



(with quark sources at infinities) these constants can be calculated as a function of the 
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Ginzburg-Landau parameter A/e 2 , see Tab. 1, where for convenience we have introduced 
the dimensionless quantities S' c = S c /M 2 and = S<p/M 2 . The numerical solution of the 
equations of motion shows that both S c and exist, are real and positive |23|. Expanding 
Eq. ( OBI ) for small x± and keeping only the leading terms, we get an equation for the 



function /: 



^/(*ll) = ^/(*ll)-^(*ll)- (4-18) 
A solution of this equation is: 

'<*> = i A '° (t) ■ (4 - 19) 

where £ = c . We remember that K (\x\/£) ~ —7 + log2 — log(|x|/£) in the short-range 



7T *L 

region (\x\ — > 0) and K (\x\/£) ~ J — - — -e~^ x ^ 1 in the long-range region (|x| — > 00). 



Since a solution exists our technical assumptions ( |4. and (4.15) are self-consistent. 
Putting Eq. ( fl~I3Q into Eq. QOTD we obtain: 



^4i4(x|| - = -9 2 /C p (a;|| - 2/,,) + ^/(x„ - yy). (4.20) 

The long-distance exponential fall off and the weakly singular (~ log(|x|)) short range be- 
haviour of the non-perturbative contribution to £1414 in Eq. (f4.20|) is compatible with the 
lattice parameterization ( |1.4|) . This fact provides an extremely interesting consistency check 
to the validity of the duality assumption (|4.1|) . Moreover this suggests the identification of 
the correlation length T g , associated with the long-range behaviour of the QCD non-local 
condensate with the dual quantity i (see Eq. ( |4.19| )). Notice that at variance with respect 



to the London limit result, here the correlation length is not simply given by the mass M of 
the dual gluon. 

Due to the almost regular short range behaviour of the non-perturbative part of Eq. 
( [4.20j ) the static potential can be calculated exactly without the use of an ultraviolet cut-off 
(at variance with respect to the London limit case, see Eq. (|3.10| )), and it is given by 
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1 „ , „ c l , C K 



V (R) = lim-5(C M , 

1 — >oo ± 



rR r+oo \ 

S c / dx 1 2{R - xi) / dxi-—-;K 

JO J~oo 2lTZ z 



Xl + xl \ g 2 1 



+ Higgs contributions 



= R £ St+ (e-«>- l )l. S .t-£} i + Rn ,4.21) 
-^Rfs c + Ra R . (4.22) 

For the simple case of = 4>(x±) the Higgs contribution to the static potential turns out 

to be given by a linear term with string tension oh-Q Taking explicitly into account this 

contribution, the total string tension is a = — S c + an = (j) 2 ) (27iS' c + a' u ) where a' K = a^/M 2 . 

2n 

For some values of a' n see Tab |. In particular, for a superconductor on the border (A/e 2 = 

1/2) from Tab. | we have Vo(R) = 7t0q [R + £ (e~ R ^ e — 1 ) ) — . In order to compare this 

^ v ' ' Air R 

potential with the heavy quark static potential we have to multiply it by the colour factor 

4 

4/3. For a typical value of O - 210 MeV we get a = -7T0 2 , ~ (430 MeV) 2 . In Fig. | we 

3 

compare the static potential of Eq. ( |4.21|) for a superconductor on the border between type 
I and type II for some typical values of the parameters with the lattice fit of Ref. ]30| . 

One of the most interesting point is to relate the dimensional parameters F 2 and T g , the 
gluon condensate and the correlation length of QCD, to the dimensional parameters 0o and 
£, the Higgs condensate and this characteristic length in the dual Abelian Higgs model. Our 
derivation identifies I with the correlation length T g and eventually explains the existence 
of a finite correlation length in terms of an underlying dual Meissner effect that gives a 



4 The comparison between Eq. ( 4.21 ) and Eq. (|2.3|) suggests to identify D(x\\) with f(x\\)S c /Tr. 
The same string tension (for what concerns the non-Higgs part) would, then, be obtained by 



using equation (2.4). We see, therefore, that the string tension is always emerging in the limit of 
large interquark distances and via an integral on a function depending on the correlation length. 
Therefore our calculation confirms the existence of the non-local condensate and traces their origin 
back to a dual Meissner effect. 



mass to the dual field. In the dual theory J2£| using trace anomaly it is possible to relate 
the Higgs condensate to the gluon condensate, F 2 ~ A0q. Using the above value of 0o and 
A/e 2 = 1/2, one obtains for the gluon condensate the value found by |[| F 2 ~ 0.013 Gev 4 . 
This is how originally it was shown in DQCD that the QCD vacuum is compatible with a 
dual superconductor on the border between type I and II JljJ . Finally, we notice that in pure 
gluodynamics the lowest dimensional gauge and Lorentz-invariant operator has dimension 
four its vacuum expectation value is the gluon condensate. In the dual model we have a 
relevant condensate, the Higgs condensate 0q, of dimension two. This could yield some 



interesting consequences in renormalon physics [31]. 

Via numerical solution of the coupled equations for C M and <p and subsequent numerical 
interpolation, it is also possible to calculate the form of all semirelativistic qq potentials 
|l^,28]. In principle, we could obtain an analytic solution also for the spin dependent and 



velocity dependent potentials following the same line of this paper. However, to this aim 
the calculation of different components of the tensor Q a ^\ p is necessary and some technical 
difficulties arise due to the fact that the simple assumption ( [4.13 ) is no longer valid. In 
the present situation we can try to gain some indications from the London limit result. 
Although, as we have seen, this is not the right limit in which to calculate the potentials, the 
qualitative long-range behaviour for the field-strength correlator is reasonable. In fact, in 
that limiting case it is possible to calculate the whole tensor G al \ p unambiguously in terms 
of some functions D and Di (see Eq. (|3.5|) ). Once we accept that in the presence of the 
quarks the short range behaviour of the Higgs field would regularize these functions on the 



flux-tube string, using the formulas of we can express all the heavy-quark potentials 
in terms of integrals over these functions D and D\. Since these functions are reasonably 
compatible with the lattice fit ( |1.4j ) this would explain the striking similarities in the long- 
distance behaviour of the potentials obtained in DQCD and in the Gaussian approximation 
of QCD 0. 

As a final remark, we notice that the flux tube structure between two heavy quarks has 
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been obtained in DQCD [TJJ as well as in the Gaussian approximation of QCD [I(J and the 
result compare very favourably in both cases with the lattice calculation ||32|| . The profile of 
the longitudinal electric field, i. e. along the string between the quarks, as a function of the 
transversal distance from the string is controlled by the penetration length in one case and 
by the correlation length in the other. 



V. CONCLUSION 

Under the assumption that the infrared behaviour of QCD is described by an effective 
Abelian Higgs model we have related the non-perturbative behaviour of the gauge invariant 
two-point field strength correlator (g 2 F fJiU (x) U(x, y)F\ p (y)U(y, x)) in QCD with the dual field 
propagator in the Abelian Higgs model of infrared QCD. In this way the origin of the nonlocal 
gluon condensate is traced back to an underlying Meissner effect and the phenomenological 
relevance of the Gaussian approximation on the Wilson loop is understood as following 
from the classical approximation in the dual theory of long distance QCD. In particular 
the correlation length T g of QCD, which we know from direct lattice measurements, can 
be expressed completely in terms of the dual theory parameters (£). As a further check we 
have calculated analytically the static potential and the string tension which are quantities 
directly related to phenomenology. It turns out that the string tension is given by an integral 
over a function of the correlation length which can be identified with the non-local gluon 
condensate. There is no cutoff introduced in this calculation since it is not performed in 
the London limit. We have shown that this limit is quite unphysical in the presence of 
sources and is valid only in the case of large distance from the chromoelectric string (which 
is different from large qq distances). Finally, these results shed some light also on the fact 
that the heavy quark potentials turn out to be equivalent in the SVM and in DQCD at large 
qq distances. 
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APPENDIX: 



In this Appendix we study the equation of motion ( |4.5|) in the presence of two static 
sources g and — g evolving from time — T/2 to time T/2 in the positions R/2 and —R/2 of 
the X\ axes respectively. Therefore x± = (£2,^3)- Under these conditions the Dirac string is 
given by: 

G%(x) = ge u , u 5 2 (x ± )(6(x 4 + T/2) - d(x, - T/2)){9{x 1 + R/2) - 9(x 1 - R/2)) (Al) 

Defining Ca(x±) = —— [ c1 2 x\\Ca(x), Eq. fl4.5|) can be written as: 
RT J 

{d\8 A B - d A d B - e 2 4> 2 (x±)5AB)C B (x±) = -ge B Aud B 8 2 (x±). 

It is convenient to split the field into the sum of two parts, Ca = C\ + C^ p , satisfying the 
equations: 

dlC p A (x ± ) = -ge BAU d B 5 2 (x ± ), (A2) 
(d 2 ± 6 AB - d A d B - e 2 <j) 2 (x ± )6 AB )C B p (x ± ) = e 2 <j) 2 (x ± )C p A (x±). (A3) 

The solution of Eq. ( |A"2^ ) is 

C A {x±) = 2 = » ( A4 ) 

e Xj_ e x± 

where we have used the Dirac quantization condition, g = 2ir/e, and 6 = (—x 3 /x±,x 2 /x±) 
is the angular unit vector in the transverse plane. Substituting (|A4|) in ( |A3| ) and defining 
C^ p (xj_) = e B Ai4,C np (x±)x B /x_i (or, which is the same, C np (x±) = C ap (x±)9, where C np is 
a vector in the transverse plane), we obtain 

(d 2 ± - e 2 <p 2 (x ± )) (^C» p (x ± )) = -e0 2 (x ± )^, (A5) 
v ' \x± J x\ 



or 



d ( 1 d (x±C ap (x ± ))) =e 2 (c^(x ± )-^j<p 2 (x ± ). (A6) 



dx± \x±dx± J \ ex±. 
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We can solve equation (HE) for small values of x±, assuming = - - + • • •, 

e 

obtaining: 

C np (* ± ) = ^x ± - 5 X 3 + . . . , (A7) 
e 8e 

where S r and S 1 ^ are some constants. 
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TABLES 



Type of Superconductor 


A/e 2 


S> 




< 


I 


1/32 


0.1125 


0.2516 


1.142 


between I and II 


1/2 


0.25 


0.6 


vr/2 


II 


2 


0.38 


1.017 


1.82 


II 


8 


0.568 


1.823 


2.06 


II 


16 


0.685 


2.49 


2.16 



TABLE I. Some values of the dimensionless quantities S' c , S'^ and a' u as a function of the 
Ginzburg-Landau parameter A/e 2 , obtained by solving the static equations of motions with quark 
sources at infinities. 
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FIGURES 



V(R) (GeV) 




R (fm) 



FIG. 1. The static potential of Eq. (4-21) for a superconductor on the border between type I 



and type II with 
Ref [30]. 



4 g 2 

210 MeV, £ = 0.22 fm and = 0.32 in comparison with the lattice fit of 

3 4-7T 



28 



